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Abstract 

We consider n particles 0 < X\ < Xi < • • • < x n < +oo, distributed according to a 
probability measure of the form 

n 

y~ n ^ n - x i) n x ? e ~ xj «> -l & > o, 

n 1 <i<j<n l<i<j<n j =1 

where is the normalization constant. This distribution arises in the context of modeling 
disordered conductors in the metallic regime, and can also be realized as the distribution 
for squared singular values of certain triangular random matrices. We give a double contour 
integral formula for the correlation kernel, which allows us to establish universality for the 
local statistics of the particles, namely, the bulk universality and the soft edge universality 
via the sine kernel and the Airy kernel, respectively. In particular, our analysis also leads 
to new double contour integral representations of scaling limits at the origin (hard edge), 
which are equivalent to those found in the classical work of Borodin. We conclude this 
paper by relating the correlation kernels to those appearing in recent studies of products of 
M Ginibre matrices for the special cases 9 = M g N. 


1 Introduction and statement of the main results 

1.1 Biorthogonal Laguerre ensembles 

The biorthogonal Laguerre ensembles refer to n particles x\ < ■ ■ ■ < x n distributed over the 
positive real axis, following a probability density function of the form 

1 n 

—A(xi,...,x n )A (xl,...,x%)JJ_±fe- x *, a > —1, 6 > 0, (1.1) 

n j =i 

where 

« n 

Z n = Z n (a,0)= / A(xi,... ,x n )A(xi,... ,x e n ) TT z“e -3 b dxj 

J[0,oo) n fJl 

is the normalization constant, and 

A(A 1 ,...,A n ) = [I (L - A *) 

1 -<i<j<n 


is the standard Vandermonde determinant. 
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Densities of the form CEO were first introduced by Muttalib [38], where he pointed out 
that, due to the appearance of two body interaction term A(aq,... ,x n )A(x 9 , ... ,x e n ), these 
ensembles provide more effective description of disordered conductors in the metallic regime 
than the classical random matrix theory. A more concrete physical example that leads to m 
(with 6 = 2) can be found in [36], where the authors proposed a random matrix model for 
disordered bosons. These ensembles are further studied by Borodin [TO] under a more general 
framework, namely, biorthogonal ensembles. It is also worthwhile to mention the work of 
Cheliotis [12], where the author constructed certain triangular random matrices in terms of a 
Wishart matrix whose squared singular values are distributed according to CEO; see also [22], 
Note that when 6 = 1, m reduces to the well-known Wishart-Laguerre unitary ensemble and 
plays a fundamental role in random matrix theory; cf. mm- 

A nice property of (11.11) is that, as proved in [38], they form the so-called determinantal 
point processes mm- This means there exits a correlation kernel K^' 9 \x,y) such that the 
joint probability density functions (11.11) can be rewritten as the following determinantal forms 

Idet 

The kernel K^' 9 \x,y) has a representation in terms of the so-called biorthogonal polynomials 
(cf. [20] for a definition). Let 

p ( '°‘ , 6 \x) = KjX j + ..., q£' e \x) = x k + ..., Kj > 0, (1.2) 

be two sequences of polynomials depending on the parameters a and 6 , of degree j and k 
respectively, and they satisfy the orthogonality conditions 

roo 

J pf ,e \x)q^ ,e \x 6 )x a e~ x dx = 5 jjk , j,k = 0,1,2,.... (1.3) 

Note that the polynomial is normalized to be monic. We then have 

n—1 

K^' d \x,y) = {x)qf' 9) {y 6 )x a e ~ x . (1.4) 

l=o 

The families {pj, j = 0,1,...} and {qk = 0,1,...}, which are called Laguerre biorthogonal 
polynomials, exist uniquely, since the associated bimoment matrix is nonsingular; see m 
and (12.61) below. The studies of these polynomials (with 6 = 2) might be traced back to [36] 
during the investigations of penetration and diffusion of x-rays through matter. Later, intensive 
studies have been conducted on the case 0 G N = {1,2,...} in [IT} [23] 23..30}03} 03}|37], where 
the general properties including explicit formulas, recurrence relations, generating functions, 
Rodrigues’s formulas etc. are derived. 

As determinantal point processes, a fundamental issue of the study is to establish the large n 
limit of the correlation kernel (11.41) in both macroscopic and microscopic regimes. By expressing 
K^' 9 \x, y) as a finite series expansion in terms of x ke y r , k, r = 0,1,..., n — 1, it was shown by 
Borodin m Theorem 4.2] that 

Kn’ 9) (^ 70 ,^ 70 ) “ {-l) k x a+k (-1 ) l y dl 6 

n 1 / 0 ~ ^ klT( a +l+ k ) lT(a + 1 + 61) a + 1 + k + 61 

k,l =0 V V / 

= 6 x a j Ja+i i(ux)J„ + i n((uv) 9 )u a du, (1.5) 

Jo e ’ e 
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( 1 . 6 ) 


where J a & is Wright’s generalization of the Bessel function m given by 

j = y- (~ X Y _. 

j ^j!r(a + 6j)’ 

see also [36] for the special case 9 = 2, a € N U {0}. These non-symmetric hard edge scaling 
limits generalize the classical Bessel kernels [20l [50 i (corresponding to 9 = 1), and possess 
some nice symmetry properties. Moreover, they also appear in the studies of large n limits 
of correlation kernels for biorthogonal Jacobi and biorthogonal Hermite ensembles [10] . When 
fl = M £ N or 1/0 = M, the limiting kernels coincide with the hard edge scaling limits of 
specified parameters arising from products of M Ginibre matrices [ 33] . as shown in [321 . 

The macroscopic behavior of the particles as n -> oo has recently been investigated in 
m , where the expressions for the associated equilibrium measures are given for quite general 
potentials and 6 > 1. According to m , as n —>• oo, the (rescaled) particles in (11.111 are 
distributed over a finite interval [0, (1 + 0) 1 + 1 / e ] ) with the density function given by 

fe{x) = i^— i ( I +(x)- I -(x)), x € (0, (1 + 9) 1+1/e ). (1.7) 

Here, I±{x) (with Irn (I + (x)) > 0) stand for two complex conjugate solutions of the equation 

/ 7 I 1\V 0 

J(z) = 9(z + l) ( —-J =x, (O,(l + 0) 1+1 / e ). 


Moreover, by m Remark 1.9], the density blows up with a rate x 1 A 1 + 61 ) near the origin (hard 
edge), while vanishes as a square root near (1 + 0) 1+1 / e (soft edge). This phenomenon in partic¬ 
ular suggests non-trivial hard edge scaling limits (as shown in (|1.5p ). as well as the expectation 
that the classical bulk and soft edge universality m (via the sine kernel and Airy kernel, re¬ 
spectively) should hold in the bulk and the right edge as in the case of 9 = 1. More explicit 

-j j n 

description is revealed later in [21]. After changing variables Xi —> 9x t , the (rescaled) particles 
are then distributed over [0, (1 + 9 ) l+e and the limiting mean distribution is recognized as 
the Fuss-Catalan distribution [3 El HQ]. Its fc-th moment is given by the Fuss-Catalan number 


1 

(1 + 0)A; + 1 


^(1 + 9)k + k^ 


k = 0,1,2,.... 


( 1 . 8 ) 


The density function of Fuss-Catalan distribution can be written down explicitly in several ways; 
cf. [|2] in terms of Meijer G-functions (see e.g. 0 EH SI] and the Appendix below for a brief 
introduction) or [34] in terms of multivariate integrals. The simplest form of the representation 
for general 9 might follow from the following parametrization of the argument 


(sin((l + 9)ip)) 1+e 7T 

sin^sin^^)) 61 ’ 1 + 0’ 


(1.9) 


It is readily seen that this parametrization is a strictly decreasing function of tp, thus gives a 
one-to-one mapping from (0,7r/(l + 9)) to (0, (1 + 9) l+e /9 e ). The density function in terms of 
<p is then given by 


P(<P) 


1 sin((l + Q)tp) . 1 (sin(/j) 2 (sin(0( / ?)) e 1 

7 tx sin (9<p) ^ 7r (sin((l + 9)ip)) e ’ 


0 < cp < 


7r 

T+9' 


( 1 . 10 ) 


From m and (11.101) . one can check directly that p blows up with a rate x near the 

origin, and vanishes as a square root near (1 + 9) 1+e /6 e , which is compatible with the changes 
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of variables. We finally note that the other description of macroscopic behavior with the notion 
of a DT-element [16] can be found in j l2] . 

The main aim of this paper to establish local universality for biorthogonal Laguerre en¬ 
sembles m- Due to lack of a simple Christoffel-Darboux formula for Laguerre biorthogonal 
polynomials, we have to adapt an approach that is different from the conventional one. The 
main issue here is an explicit integral representation of K^’ 8 \ Our main results are stated in 
the next section. 


1.2 Statement of the main results 

Our first result is stated as follows: 

Theorem 1.1 (Double contour integral representation of K^fi' 8 ^). With defined in (jl.4jl . 

we have 


Kt’ e \x,y) = 


rc+ioo 


(27 Ti) 


ds ® d t 


T(s + l)T(o: + 1 + 0s) T(t — n + 1) x 9s 1 y 9t 

r(t + i)r(o: +1 + ot) r(s — n +1) s — t 


, ( 1 . 11 ) 


for x, y > 0, where 


c = 


max{0,1 - 3+1} - 1 


< 0 , 


( 1 . 12 ) 


and T, is a closed contour going around 0,1,..., n — 1 in the positive direction and Ret > c for 
t € X. 


We highlight that this contour integral representation bears a resemblance to those appearing 
recently in the studies of products of random matrices ng m \m E3], where the integrands 
of double contour integral representations for the correlation kernels again consist of ratios of 
gamma functions. When 6 € N, Kjf 31 ’ 8 ^ is indeed related to certain correlation kernels arising 
from products of Ginibre matrices; see Section [3] below. We also note that, in the context of 
products of random matrices, the correlation kernels can be written as integrals involving Meijer 
G-functions, for biorthogonal Laguerre ensembles, however, it does not seem to be the case for 
general parameters a and 6. 

An immediate consequence of the above theorem is the following new representations of 
hard edge scaling limits. 

Corollary 1.2 (Hard edge scaling limits of K^fi’ 8 ^). With a > —1, 8 > 1 being fixed, we have 


K, 


(a,e) 


lim ■ 

n—too 


x y 
,i/e j n i/e 


n 


1/0 


= K^ e \x,y), 


uniformly for x, y in compact subsets of the positive real axis, where 

Q /*C +200 


r\ pc~\~ioo r 

K {a ’ 8) (x,y) = ——/ d sf 

(2mY J c _ ioo ,/ s 


df 


T(s + l)T(a + 1 + 8s) sin7rs x ds l y et 
T(t + l)T(a + 1 + 8t) sin7rt s — t 


(1.13) 


(1.14) 


and where c is given in (11.121) . E is a contour starting from +oo in the upper half plane and 
returning to +oo in the lower half plane which encircles the positive real axis and Re t > c for 
t € H. Alternatively, by setting 


p{oL,e) ^ 


j_ r +io ° x s dc 

2W™ T((l — s)/6) 


k > —a , 


(1.15) 
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and 


(1.16) 


Q 


( a > fl >(s) = — f 

K J 2 m L 


nm 


x f d t, 


,ry r(Q! +1 - t) 

where 7 is a loop starting from —00 in the lower half plane and returning to —00 in the upper 
half plane which encircles the negative real axis, we have 


K^ e \x,y) = f 1 p^ e \ux)q^ e) {uy)du. 

Jo 


(1.17) 


In Corollary 11,21 we require 8 > 1 to make sure the integral is convergent. Note that when 
0 = 1, we have (see jjlj formula 10.9.23]) 

p^Xux) = ( ux) a/2 J a (2y/ux ), q^Xuy) = (uy)~ a/2 J a (2y/uy), 

where J a denotes the Bessel function of the first kind of order a. It then follows from (11.171) 
that 


K^ 1 \x,y) = 


a/2 ,1 


J a (2^/ux)J a (2y/uy) du 


= 4 I — 


a/2 


K?*(4x,4y), 


where 




Ja(Vx)^yJ'JVy) - VzJL('/x)Ja(Vii) 


a > —1, 


2(x - y) 

is the Bessel kernel of order a that appears as the scaling limit of the Laguerre unitary ensembles 
at the hard edge [20,150], as expected. Furthermore, a comparison of (11.51) and (11.131) (11.141) 
gives us the following identity 


8x a [ J a +1 1 (ux)J n +i v((uy) e )u a du 

Jo e ’ 6 

pc-\-ioo 


q rc-\-i 00 p 

Fi) 2 Jc-ioo d ' S Jz 


(27TZ 


dt 


r(s + l)r(d! + 1 + 9s) sin7rs x 9s 1 y et 
r(t+l)r(a + l + 0t) sinvrt s — t 


> 1. (1.18) 


For a direct proof of the above formula; see Remark 13.21 below. 

We believe that the new integral representations (11.141) and (11.171) for A"(a,0) w m a j so facili¬ 
tate further investigations of relevant quantities, say, the differential equations for the associated 
Fredholm determinants, as done in m-m The studies of these aspects will be the topics of 
future research. 

By performing an asymptotic analysis for the double contour integral representation (11.111) . 
we are able to confirm the bulk and soft edge universality for biorthogonal Laguerre ensembles, 
which are left open in [13] . The relevant results are stated as follows. 

Theorem 1.3 (Bulk and soft edge universality). For xq £ (0, (1 + 8) 1+d /8 9 ), which is parame¬ 
terized through OK by ip = p(xq) € (0,7r/(l + 8)), we have, with a, 8 being fixed, 


p —nr) cotip 1 

lim ;--- r 

n —>00 e 7r, f co t ( ^ , v l—o 

pmx 0 


Ki a ’ e) \ne(x 0 + 


z 


np(ip) 


,n0 ( x 0 + 


np(p) 


= K S in (£, 77 ), (1-19) 
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Figure 1: The contours 7 l and 7 # in the definition of Airy kernel. 


uniformly for £ and rj in any compact subset of M, where p((p) is defined in (11.1011 and 

sin 7 r(x — y) 


K s i n (x, y) .— 


tt(x - y) 


( 1 . 20 ) 


is the normalized sine kernel. 
For the soft edge, we have 


lim + 1 


n ^°° e -2~^(l+0)§|n3 23 


2 

n 3 


n I . C *V 
nt) [ x* -\ -j- 

ns 


uniformly for £ and 7 in any compact subset of K, where 


x* = 


(1 + 0 ) 


1+0 


c* — 


(1 + 0 ) 


+0 


2S0 0 - 1 


= K M (S,V) (1-21) 


( 1 . 22 ) 


and 


K Ai (x,y) := 


Ai(x) Ai'(y) - Ai'(x) Ai(y) 


d/i / dA- 


3 x i l 1 


® - V (2m) 2 J yL ^' e ^-y\fi-X 


(1.23) 


is f/ie Airy kernel. In (11.2311 . 7 ^ and 7 l are symmetric with respect to the imaginary axis, and 
7 ij is a contour in the right-half plane going from e _7r//3 * • 00 to e n / 3 * • 00 ; see Fiynre [7] /or an 
illustration. 


In the special case 8 = 2, a € N U {0}, the bulk universality is first proved in [ 36] . 

Remark 1.1. The result of soft edge universality (11.211) also implies that the limiting distri¬ 
bution of the largest particle in biorthogonal Laguerre ensembles, after proper scaling, converges 
to the well-known Tracy-Widom distribution [ 6 j Theorem 3.1.5]. 
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1.3 Organization of the rest of the paper 

The rest of this paper is organized as follows. Our main results are proved in Section [2j The 
proofs of Theorem 11.11 and Corollary 11.21 are given in Sections 12.21 and 12.31 respectively, which 
rely on two propositions concerning the contour integral representations of and in 

Section \2 .11 These formulas might be viewed as extensions of the intensively studied 6 € N case, 
and we give direct proofs here. The nice structures of these formulas then allow us to simplify 
COD into a closed integral form as well as to obtain the hard edge scaling limits, following 
the idea in recent work of the author with Kuijlaars [33j. The bulk and soft edge universality 
stated in Theorem 11.31 is proved in Section [2.41 We will perform a steepest descent analysis of 
the double contour integral (11.111) . whose integrand constitutes products and ratios of gamma 
functions with large arguments. It comes out that the strategy developed by Liu, Wang and 
the author in [35] (see also m) works well in the present case. Roughly speaking, the strategy 
is to approximate the logarithmics of the gamma functions by elementary functions for n large, 
which play the role of phase functions. There will be two complex conjugate saddle points in the 
bulk regime, corresponding to the sine kernel, while in the edge regime these two saddle points 
coalesce into a single one, which leads to the Airy kernel. A crucial feature of the analysis is to 
construct suitable contours of integration with the aid of the parametrization (11.91) . Since the 
asymptotic analysis is carried out in a manner similar to that performed in [35], emphasis will 
be placed on key steps and demonstration of basic ideas in the proof of Theorem 11.31 but refer 
to [35] for some technical issues. 

We finally focus on the cases when 6 = M € N, and relate to correlation kernels of 

specified parameters arising from products of M Ginibre matrices. Some remarks are made in 
accordance with this relation to conclude this paper. For convenience of the reader, we include 
a short introduction to the Meijer G-function in the Appendix. 


2 Proofs of the main results 

2.1 Contour integral representations of p^ ,6 \x) and q[ a ' e \x ) 
Proposition 2.1. We have for x > 0, 




(-x)i 


3=0 

T(a + l + k6)k\ r 


jj T(a + 1 +jd) 


—r(a + l + k6) 


T(t — k)x t 


( 2 . 1 ) 


2 ni T(t + l)r(a; + 1 + Ot ) 

where £ is a closed contour that encircles 0,1 ,..., k once in the positive direction. 

Proof. The first identity in (12.11) follows from the determinantal expressions for the polynomials 
qk l - d \ gy getting the bimoments 

r*+oo 

/o 


rn j,k ~ 


r+oo 

/ x a+j+ek e~ x dx = T(a + j+ kd+ 1), j,ke NU{0}, 

Jo 


( 2 . 2 ) 
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we define 


D n — det(rrij,k)j,k=o,...,r 

( r(a + l) 


= det 


r(a + 2) 


r(a + 1 + 0) 

r(a + 2 + 0) 


T(a + 1 + nQ) ^ 
r(a + 2 + nO) 


(2.3) 


yr(a + 1 + n) r(a + 1 + n + 0) ••• T(a + 1 + n(0 + 1 ))/ 

From the general theory of biorthogonal polynomials (cf. j!5l Proposition 2]), it follows that 

\ 


«£"*(*) = 



^0,0 

m 0t i 

TYl0,k 

-det 




D k -1 

mk- i,o 

mk- 1,1 • 

TTlk—l,k 


V 1 

X 

x k 


/ 

/r(a + l) r(a + l + 0) ••• r(a + 1 + k6) \ 


Du 


■ det 


fc-i 


T(a + k) r(a + k + 6) 

y 1 x 


r(cr + k{ 1 + 0)) 

„k 


, k > 1 (2.4) 


with q^ a,e \x) = 1. With the aid of functional relation 

r(z + 1) = zT(z), 

an easy Gauss elimination process gives us 


D n = ]JkW k r(a + l + k9). 


k=0 


(2.5) 


( 2 . 6 ) 


Similarly, by expanding the matrix in (|2.4D along the last row and evaluating the associated 
minors, it follows 


i a ’ e \x) 



( -x) j 

r(a + 1 + j9) 


r(ct! + 1 + k9), 


(2.7) 


see also [301 for a proof of (12.71) by checking the orthogonality directly if 0 = M. 

To show the second identity in (12.11) . we note that integrand in the right-hand side of (ED 
is meromorphic on C with simple poles at 0,1,... , k (the poles of the numerator at the negative 
integers are canceled by the poles of the factor T(t + 1) in the denominator). Hence, by the 
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residue theorem and a straightforward calculation, we obtain 


T(a + 1 + k6)kl f 


T(t — k)x f 


27 vi J-£ T(t + l)r(a; + 1 + Ot) 


d t 


r(a + 1 + k6)k\ Res t=j 

3=0 
k 

T(a + l + k9)klY^ 


T(t - k ) 


v(t + i)r(a +1 + et) 


3=0 


(k-j)\jW(a + l+je) 


H)"E 


{-x) j 


—T(a + l + kd). 


j=0 JJ r ( a + 1 + J e ) 

This completes the proof of Proposition 12.11 

Proposition 2.2. For p^’ 9 \ we have the following Mellin-Barnes integral representation 

i r c+io ° r(f+ 1 -|) 


Ot —X / \ 

x e p k {x) = 


27rir(a + 1 + kd)k\ „/ c _ ioo T (| + 1 — | — A:) 
where c > max{—a, 1 — 0} and x > 0. 


T(a + s)x S ds, 


( 2 . 8 ) 

□ 

(2.9) 


Proof. Note that all the poles of the integrand lie on the left of the line Re z = c, it is readily 
seen that the integral formula in the right-hand side of (12.91) is well-defined. On account of 
the uniqueness of biorthogonal functions, our strategy is to check the integral representation 
satisfies 


the orthogonality conditions 

i roo pc+ioo t (s i -1 _ 1\ 

J_ / x jo / 1 U + 1 e) 

2TriF(a + 1 + k9)k\ J 0 j c _ 

for j = 0, 1, ... , fc; 


r (f +1 - i - fc) 


T(a + s)x s ds dx = 6j k , (2.10) 


the integral ^ fc-i^ i 9 {•) b(a + s ) x S belongs to the linear span of x a e x , 




To show (12.101) . we make use of the inversion formula for the Mellin transform and obtain 

1 roo rc+ioo r /s i 1 _ 1 \ 

i / ,-A / 1 U "r 1 e) 


roc re 

-d *°L 


2TriY(a + 1 + kO)k\ 

r (I + 1 -1) 


ru, + i - kn)kw f j + i - i - k) 
(j + 1 — fc)*r(l + a + j6) 


r + i-J-i) 

T(a + s) 


T(a + s)x S dsdx 


s=jO +1 


T(a + 1 + k9)k\ 


= 5. 


i+- 


( 2 . 11 ) 


To check the second statement, recall the Pochhammer symbol (a) k = = a ( a + 

1) • • • (a + k — 1), it is readily seen that 


r(f + i-a) 

r(| + i-i-/c) 


= t +1 - t - fc = ^ + i- x- fc---- + i- -- i 
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is a polynomials of degree k in s, the integral is then a linear combination of weights Wj a \x), 
j = 0 ,k, where 


w 


i rc+ioo 

“ (x) = — / s- J r(a + s) 

Z'Kl Jc—ioo 


x S ds. 


( 2 . 12 ) 


Thus, it suffices to check Wj°\x) belongs to the linear span of x a e X ,x a+l e x ,..., x a+ ^e x . 
We now expand the monomial s J in terms of the basis (a + s)i, l = 0,..., j, i.e., 


s J = 


^ ~2ai(a + s)i = Y2 ai 


T (<y -)- l -)- s) 
T(a + s) 


1=0 1=0 

for some constants ai with aj = 1. Inserting the above formula into (12.121) . it follows that 


(a) ( \ 

Wj '(x) = 


1 J rc+ioo J 

^Ti^2 ai j c _ ioo r(a + l + s)x~ s ds = J2 


= > a lX a+l e~ x , 


(2.13) 


1=0 c-ioo 

as desired, where we have made use of the fact that 

^ rc+ioo 
2n i 


rc+ioo 

- T(v + s) 

L J C —200 


x s ds = aTe 


1=0 


v > —1; 


see (1A.4I) below. 

This completes the proof of Proposition 12.21 
2.2 Proof of Theorem 11.11 


□ 


With a change of variable s —> 9s + 1 — 9 in (12.91) and contour deformation, we rewrite 

cy —t / \ 

x e p k ( x) as 


9x e 


-l 


rc+ioo 




-T(0s + l-6 + a) X - es d s, 


2niT(a + 1 + k6)k\ J c -ioo T (s — k) 
where c > max{0,1 — ^p-}. This, together with (11.41) and (12.11) . implies that 


(2.14) 


f) T 0- 1 rc+ioo r 

Kn’ e) (x,y) = ——^ / ds A dt 

(2m) 2 J c _ ioo 


f T(s)r(0s + 1 — 9 + a) ^ F(f — k ) 0 S 0 t 

/ -j t'/ „ 7„\ ^ y 

k =0 


(2.15) 


T(t + l)T(a + 1 + Ot) ^ r(a - k) 

We now follow the idea in [33l. From the functional equation (12.51) . one can easily check that 

, N r (t-k) r(t-k) r(t-k + 1 ) 

<- t - 1) r|r^ = r(.-t-i) - r ( .-t) ■ 

which means that there is a telescoping sum 


71—1 




k =0 


T(t — k ) T(t-n + l) r(t + 1 ) 
T(s-k) = T(s-n) r(s) ' 


(2.16) 


To make sure that s-t - 1 / 0 when s € c + iM. and f G S, we make the following settings. 
Note that max{0,1 — ^p-} < 1 for a > — 1 and 9 > 0, we take 

1 + max{0,1 — ^Tp-} 
c =-^- e — < 1 
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and let X go around 0,1,... ,n — 1 but with Ret > c — 1 for t € X. Then we insert ([2. 161) into 
(12.151) and get 


M (x v) [ c+io ° ds I r(8)r(03 + l-6 + a)r(t-n + l) 

n 1 ,y) (27T i) 2 J c _ ioo T(t + l)T(a + l + 9t) T(s — n) s-t-1 

_ ex e ~ l f c+io ° /■ r (ds + 1 -9 + a) x- 8s y 8t 

(2m) 2 J c _ ioo S / s t T(a + 1 + Ot) s-t- 1' 

The f-integral in the second double integral vanishes due to Cauchy’s theorem, since there are 
no singularities for the integrand inside X. With a change of variable s i —> s + 1 in the first 
double integral, we obtain (11.111) 

This completes the proof of Theorem 11,11 


2.3 Proof of Corollary 11.21 

The proof now is straightforward by taking limit in (11.111) . as in [33] . Recall the reflection 
formula of the gamma function 


T(i)T(l — t) = — 


TT 


it is readily seen that 


sin irt' 

T(t —n + 1) T(n —s)sin7rs 
T(s — n + 1) T(n — t) sin7rf 

As n -+ oo, we have (cf. [3T1 formula 5.11.13]) 

r (™- ») _ t-s 


T(n - t) 


= n 


(l + 0(n x )) , 


(2.17) 

(2.18) 

(2.19) 


which can be easily verified by using Stirling’s formula for the gamma functions. By modifying 
the contour X in (11.111) from a closed contour around 0,1,... , n — 1 to a two sided unbounded 
contour starting from +oo in the upper half plane and returning to +oo in the lower half plane 
which encircles the positive real axis and Ret > c for t € X, the scaling limits (11.141) follow. 
The interchange of limit and integrals can be justified by combining elementary estimates of 
the sin and gamma functions with the dominated convergence theorem, as explained in [83]. 

To show (I1.17p . we note that 


-6>s-i te /• i 

---— =—9 / (ux)~ 8s ~ 1 (uy) 8t du, (2.20) 

s-t Jo 

and, by (12.171) . 

sin7rs T(1 + f)T(—t) 
sin7rt T(l + s)T(— s) 

Inserting the above two formulas into (11.141) . it is readily seen that 


K^ e \x,y) = - 


f (l, j 

JO \ /7rz Jc-io o 


r(a + 1 + 9s) 

r(-sj 


ux) 


—8s —1 


ds 


x 


9 

2ixi 


r(-t) 


T(a + 1 + Ot) 


(uy) 8t d t d u. 


The change of variables s i->- 9s + 1 and 1 1 -+ —Ot takes the two integrals into the two functions 
p( a ’ 8 ' 1 and q^ a ' e "> defined in (11.151) and (11.161) . respectively. The identity ([1.1711 then follows. 
This completes the proof of Theorem 11.21 
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2.4 Proof of Theorem 11.31 

We start with a scaling of the correlation kernel K^’ 6 \x, y ) —> K^'°\0x^ ,9yi>). By (11.111) . 
it then follows that 


i4 a ’ e \0xKoy*) 

_ i r ds r d ^ r(s + i)r(a + 1 + o 8 ) r(t - n + 1) g-^a-^y 
(2iri) 2 xi> Jc %: B(t + l)r(a + 1 + 6t) T(s — n + 1) s — t 


( 2 . 21 ) 


where C and £ are two contours to be specified later, depending on the choices of reference 
points. 

By setting 


F(z-a) := log 


/T(z + l)r(a + 1 + Oz) dz \ 
{ r ( z-n + i)- s “ j 


a > 0, 


( 2 . 22 ) 


where the branch cut for the logarithmic function is taken along the negative axis and we assume 
that the value of logz for z € (—oo, 0) is continued from above, we could rewrite (I2.2ip as 


ht’ e \exKe y *) 


i) 2 xe Jc Jt, 


(: 2ni) 2 xe Jc 


e F(s;x) y 
dt e F(f,y) s _ t - 


(2.23) 


We will then perform an asymptotic analysis of (|2.23l) . The basic idea is the following. It 
is clear that the function F in (12.231) plays the role of a phase function. For large 2 and proper 
scalings, F can be approximated by a more elementary function F (see (12.291) below) with the 
help of the Stirling’s formula for gamma function. There will be two complex conjugate saddle 
points w± (see (12.321) below) of F in general. In the proof of bulk universality, the two contours 
are deformed so that one of them will meet the pair of saddle points. It comes out that the 
main contribution to the integral does not come from the saddle points alone, but from the 
vertical line segment connecting the two points. In the proof of soft edge universality, the two 
saddle points coalesce into a real one. The phase function then behaves like a cubic polynomial 
around the saddle point (see (12.471) below), which justifies the appearance of Airy kernel. 

We also note the possibilities to deform the contours in (12.231) . Firstly, it is readily seen that 
the integral contour for s can be replaced by any infinite contour C oriented from —zoo to zoo, 
as long as £ is on the right side of C. One can further deform C such that £ is on its left, and 
the resulting double contour integral remains the same. To see this, let C and C be two infinite 
contours from —zoo to zoo such that £ lies between C and C. An appeal to the residue theorem 
to the integral on CuC gives 



d t 


e F{s\x) ^ 
e F(t ]y ) s _ t 




p F{s\x) 1 

d , 

e F(t-,y) s _ f 



d t = 0. 


(2.24) 


Hence, the double contour integral does not change if C is replaced by C. We will use such kind 
of contour deformation in the proof of the soft edge universality. Similarly, one can show that if 
£ is split into two disjoint closed counterclockwise contours £ = £i U £ 2 , which jointly enclose 
poles 0,1,..., n — 1, and C is an infinite contour from —zoo to zoo such that £1 is on the left 
side of C and £2 is on the right side of C, the formula (12.231) is still valid. We will use such kind 
of contours in the proof of the bulk universality. 

We now derive the asymptotic behavior of F. Recall that the Stirling’s formula for gamma 
function [0Tj, formula 5.11.1] reads 


log r(z) 



z + ~ log(2vr) + O 



(2.25) 
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as z —>■ oo in the sector |arg£j < tt — e for some e > 0. It then follows that if |z| 
\z — n\ —>• oo, while arg 2 and &rg(z — n) are in (—n + e, ir — e), then uniformly 

F(z ; a) = F(z ; a) + ^(log 2 - log (2; - n)) + ^ log(27r) + 0(min(|z|, |z - nl)" 1 ), 
where 

F{z\a) = (1 + 9)z(\ogz — 1) — (z — n)(log(z — n) — 1) — z log a. 
Furthermore, we have 

F(nz ; n e a ) = nF(z\ a) + n logn, 

where 

F(z\ a) = (1 + 0)^(log z — 1) — (z — l)(log( 2 ; — 1) — 1) — zloga. 

Note that if 6 = M e N, we encounter the same F(z; a) and F(z; a) as in [35]. 

Since 

F z (z] x) = (1 + 9)logz — log (2 - 1) - log®, 
the saddle point of F{z\x) satisfies the equation 

z 1+e = x(z - 1). 


00 and 

(2.26) 

(2.27) 

(2.28) 

(2.29) 

(2.30) 

(2.31) 


In particular, if x = xq € (0, (1 + 9) l+d /9 e ), which is parameterized through (11.91) by tp = 
tp(x 0 ) € (0, 7 r/(l + 9)), one can find two complex conjugate solutions of (12.311) explicitly given 
by 

= sin((l + fl)¥?) ±i 
sin(0(^) 


(2.32) 


For later use, we also define a closed contour 

S - L - Sin (( 1 + 6 W M 


sin (9(j)) 


IT IT 

< 6 < 


1 + 9 


1 + 9 r 


(2.33) 


which passes through w±, intersects the real line only at 0 when (f> = +ir /(1 + 9) and at 1 + 9~ 1 
when (p = 0 . Since the integrand of (12.231) takes 0 as one of the poles, we further deform £ a 
little bit near the origin simply by setting 

£ e := {z € £ | \z\ > e} U the arc of {|z| = e} connecting £ D {|^r| = e} and through —e, (2.34) 
with counterclockwise orientation. 

With the above preparations, we are ready to prove the bulk and soft edge universality for 
K (a,e) 


Proof of (11.191) In view of (11.191) . we scale the arguments x and y in ([2.231) such that 


x = n I Xq + 


np{ip) 


y = n [x 0 + 


np(<p) 


(2.35) 


where £ and 77 are in a compact subset of R and p{ip) is given in ( 11 . 101 ) . 

The contours C and £ are chosen in the following ways. The contour C is simply taken to be 
an upward straight line passing through two scaled saddle points nw±. This line then divides 
n£ r into two parts, where r is a small parameter depending on 9. By further separating these 
two parts, we define 

£ = £cur u £ V erj (2.36) 
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Figure 2: The contours C and £ used in the proof of bulk universality. 


where £ cur is the part from n£ r , and £ ver are two vertical lines connecting ending points of 
£ cur . The distance of these two lines is taken to be 2e, with C lying in the middle of them; 
see Figure [2] for an illustration. The main issue here is that, with these choices of C and £, 
ReF(z;xo) defined in (|2.29l) attains its global maximum at z = w± for nz € C and its global 
minimum at z = w± for z € £, which can be proved rigorously with estimates as shown in [33J 
Lemma 3.1]. 

By taking the limit e —» 0, it follows 

K^ e \exK0y^ = h + I 2 , (2.37) 

where (p. v. means the Cauchy principal value) 


11 := lim-- 

e_s, ° [2fi) 2 x~s 



df 


e F(s-,x) 

e F ^y) 


1 


s — t 



gF(s,x) 


1 


s — t 


d t, 


(2.38) 


and, by interchange of integrals and the Cauchy’s theorem, 


I 2 : = lim- 

e_5, ° {2m) 2 x^ JC 


1 1 ds I 

0 JC J £, 


d t 


2 irix 


IX 6 Jn 


e F{s-,x) ^ 

e F ^y) s -1 

nw+ 


2 irix 


^ j 

IX 6 Jn 


nw-y. F(s\x) 


?F(s\y) 


ds 


^Vds = 


2 irixe log(|) 


y\nw+ 

x) 



(2.39) 


Here we note that by taking e —>• 0, the vertical line (nw-,nw + ) is enclosed by £ ver , hence the 
Cauchy’s theorem is applicable in the first step. 

With the values of x,y given in (12.3511 and w± given in (12.3211 , a straightforward calculation 
gives us 


I 2 = 


p(¥>)®o ® 


2 m(ri - £) (1 + 0 (n -1 )) 

r_i e 7rcot w S i n Tr^-jy) 


(v-i)ni + 

e pM *0 (l + O (n 


-1 


))- 


e PM *0 (l + O ( 


-1 


n 


)) 


= P(v)x 0 


g 7 T COt ip£ 7]-(^ — 77) 


+ 0 (n~ l ) 


(2.40) 
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Figure 3: The contours C and £ used in the proof of soft edge universality. 


On the other hand, one can show that, in a manner similar to the estimates in | 35l Lemma 
2.1], F{z\n e x o) attains its global maximum at z = nw± for z € C and its global minimum at 
2 = nw ~t for z € £, which leads to the fact that I\(z) = 0(n -1//2 ). This, together with (12.3711 
and (12.401) . implies (|1.19l) . 


Proof of (11.211) On account of the scalings of x,y in (11.211) . we set 

x = nd ( x * + ’ y = n e (x* + , (2.41) 

where £,ry S M, x* and c* are given in (|1.22l) . 

In this case, the two saddle points w± coalesce into a single one, i.e., 

w + = w- = zo = 1 + -. (2.42) 

U 

We now select the contours £ and C as illustrated in Figure [3] The contour £ is still a defor¬ 
mation of n£ r , while near the scaled saddle point nzo, the local part £i oc is defined by 


Sioc = jnz 0 + cin3re 2m/3 r€ l,nso | U |nzo + cinsre 2m l‘ i re l,nso | 


U 


2 

cins 2 

nzo-1- icin^r 


r £ 


V3 Vz 


where 

ci := x*/c* = 25(1 + 0)5/0- (2.44) 

The contour C is obtained by deforming a straight line. Around nzo, the local part is defined 
by 


Cioc = 1 nz 0 + Cinzre 


2 


r € 


1 , 77,30 


M 


nzo + cin^re 


Sl-P Tti/'i 


c\n3 . 2 

U < nzo H-1- icin^r 


r € 


r € 


1,7130 


V3 v/3 

IP IT 


(2.45) 


As in [351 Equation (2.69)], one can show that the main contribution to the integral ([2.231) . 
as n —>■ oo, comes from the part C\ oc x £i OC) and the remaining part of the integral is negligible. 
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When (s,t) € C\ oc x Xi oc , we can approximate F(s-,n d x ie ) and F(t,n e x *) by F given in (12.261) . 
and further by F that is defined in (12.291) . 

With zg given in (12.421) . it is readily seen that 


F z (zg, x*) — 0, F zz (zg, X*) — 0, F 22;2 ;(zo, 3?*) — 


0 3 

1 + 0 ' 


(2.46) 


Hence, 

F(zo + n~3c\u; x*) 

i i „ 2 1 - 

= F(z 0 ;x *) + F z (zg-,x*)c\urTz + -F zz (^ 0 ;x*)c)Vn _ 3 + -F 222 (z 0 ;x*)cfu 3 n" 1 + 0 (n“5 


— F( z o'i x *) + ^ + 0 (n 5 ^ • 

By changes of variables 

2 2 
s = nzo + ^ 3 ciu, t = n^o + ^ 3 ciu, 

it follows from (12.231) . (12.261) . (12.411) . (12.441) and (12.471) that 

e F(s-,x) i 


(2.47) 

(2.48) 


1 


ds 


i 


d t- 


(2m) 2 xe JC loc 7s loc e F ^s-t 


(2ni) 2 xo JCi 0 


e F(s;n M x *) ^ _)_ n 3 q ^ ^ 


ds ® df 

'Si oc e F(t\n M Xt) ^1 _|_ 


—t 


s — t 


o 2-3(l+0)3(r?-e)n^ 


Cl 


i 

nsxS 


(2vri) 2 J c dU f T 


dw 


ea 


ilI 3 -U^ 


Co J Sq e s 


\v z -vq u — V 


+ 0 [n 5 


2 § e 2 - 3 ( l +++ 7?-£+ 3 

~ 1,2 

n3 (1 + 6)o + 3 


KAi(£,v) + ° ®)) > 


(2.49) 


where So and Cg are the images of C\ oc and Si oc (see (12.451) and (12.431) 1 under the change of 
variables (12.481) . and the last equality follows from the integral representation of Airy kernel 
shown in (11.231) . 

This completes the proof of Theorem 11.31 


3 The cases when 0=MeN 

In this section, we will show a remarkable connection between K and those arising from 
products of Ginibre random matrices if 0 € N. In the limiting case, this relation has been 
established in [32]- Our result gives new insights for the relations between these two different 
determinantal point processes. In particular, it provides the other perspective to explain the 
appearance of Fuss-Catalan distribution in biorthogonal Laguerre ensembles; see Remark 13.11 
below. We start with an introduction to the correlation kernels appearing in recent investigations 
of products of Ginibre matrices. 

3.1 Correlation kernels arising from products of M Ginibre matrices 

Let Xj, j = 1,... ,M be independent complex matrices of size (n + Uj) x (n + Vj~\) with 
ug = 0 and Vj > 0. Each matrix has independent and identically distributed standard complex 
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Gaussian entries. These matrices are also known as Ginibre random matrices. We then form 
the product 

Ym = XmXm- i ■ ■ ■ X\. (3.1) 

When M = 1, Y\ = X\ defines the Wishart-Laguerre unitary ensemble and it is well-known 
that the squared singular values of Y\ form a determinantal point process with the correlation 
kernel expressed in terms of Laguerre polynomials. Recent studies show that the determinantal 
structures still hold for general M M- According to [3], the joint probability density function 
of the squared singular values is given by (see 0 formula (18)]) 


P(x i, ...,x n ) = —A(xi, ...,x n ) det [wk-i{xj)\ jk=l jl , Xj > 0, 


(3.2) 


where the function w k is a Meijer G-function 


w k (x) 


G, 


M, 0 
0 ,M 


I'M, VM- 1, ■ ■ ■ ,V2,Vl 



(3.3) 


and the normalization constant (see [3, formula (21)]) is 


n M 

Zn = n\ JJr(i + Uj). 
i=l j=0 


Note that the Meijer G-function w k (x) can be written as a Mellin-Barnes integral 


^ pc+io o M 

w k{x) =-— / T(s + ui + k)W r(s + Vj) X - s ds, A: = 0,1,..., (3.4) 

Jc —zoo j _2 

with c > 0. As a consequence of (IA.4D . it is readily seen that if M = 1, (13.21) is equivalent to 
(11.11) with 6 = 1. 

The determinantal point process (13.21) again is a biorthogonal ensemble. Hence, one can 
write the correlation kernel as 


71—1 


K(x, y ) = Y.n<.x)Ql(y), 

k=0 


(3.5) 


where v stands for the collection of parameters v\., vm and the biorthogonal functions PY 
and Q v k are defined as follows. For each k = 0,1,... ,n — 1, P k is a monic polynomial of degree 
k and Q v k can be a linear combination of wo,... ,Wk, uniquely defined by the orthogonality 

Pj (x)Q k (x) dx = S jik . (3.6) 

In particular, we have the following explicit formulas of P k and Q k in terms of Meijer G-functions 

0: 



Qfc(z) = 


nf= 0 nk+^ + i) 


r M+1,0 ( ~k 
( “T,M+l 1 ,, „ 

V^O ,V\,...,VM 


1 


2?r ^ nflo + !) Jc- 


G=o 


pc- 1 - 
J c—i 




T(s-k) 


(3.7) 


17 









and 


M 

P^x) = -l\T(n + ^ + l)G^ M+1 

j =o 

_ ( nnrr r ( n + v j + !) p 

-(-D n T(yj + 1) ,F " 


n + 1 

-VQ, —V 1 , • • • , —I'M- 1 , —VM 


—n 

1 + z/l, - - - , 1 + I'M 


X 


(3.8) 


where 


/ ftl, . . • , (Ip \ _ \ (ftl)fc ' ’ • (ftp)fc ^ 

\ A / (^l)fe • • • {bq)k k\ 


is the generalized hypergeometric function with 


(a)k 


T(a + k) 

r(a) 


a(a + 1) • • • (a + k — 1) 


(3.9) 


(3.10) 


being the Pochhammer symbol; see (IA.2D for the second equality in (13.81) . The polynomials 
P k can also be interpreted as multiple orthogonal polynomials |23] with respect to the first M 
weight functions Wj, j = 0,... ,M — 1, as shown in [33]. More properties of these polynomials 
(or in special cases) can be found in [H[ 33 3D] [52] E>3> 55]. 

With the aid of (13.71) and (13.81) . it is shown in [331 Proposition 5.1] that the correlation 
kernel admits the following double contour integral representation 




1 r 1/2+i °° f yr T(s + Vj + 1) T(t - n + 1) xV*" 1 

(27T*) 2 J-l/2-ioo Ve 7=0 P(* + V i + 1) P( s - ^ + 1) S-t ’ 


(3.11) 


where S is a closed contour going around 0,1,...,n — lin the positive direction and Re t > —1/2 
for t € X. For recent progresses in the studies of products of random matrices; see [2]. 

We point out that the kernel (13.111) (as well as the biorthogonal functions P k and Q k ) is 
well-defined as along as i'i > —1, i = 1 ,... , M, and has a random matrix interpretation if z/j 
are non-negative integers, i.e., then the particles correspond to the squared singular values of 
the matrix Ym- 


3.2 Connections between Ki oM> and Kf 


Our hnal result of this paper is stated as follows. 

Theorem 3.1 (Relating Kn i,M> to K”). Let p^ k ' e \ q k ' 6 ' > , P k and Q k be the functions defined 
through biorthogonalities (USD and flSBD , respectively. If 8 = M we have 


q { k a,M \x) = M kM Pfi 


k V M m ) ’ 


X 


M—l 


x'e-'p^x) = Ql (^g) M( ; +1)M _ r 


where the parameter D is given by an arithmetic sequence 

QL 7 

Dj = - 1-1, 7 = 1,..., M. 

3 M M ' J ’ ’ 

As a consequence, we have 

/ y M x M \ 


M -1 


K^n'tfx y) = — _ K° ( —_-_ ] 

n j^m-i n I M M ’ M M J 


(3.12) 


(3.13) 


(3.14) 


where 


and K” are two correlation kernels defined in OH and ra, respectively. 
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Proof. Suppose now the parameters in P£ are given by (13.131) . we see from (13.81) that 

Hi?) 


k l r + 1 ) 


'i'll 


E 


(-fc)t 


x 


j=i + 1) i=0 (1 + ^l)i‘ ‘' (1 + vm)% *! 


M 


(-i) fc n r ( fc + 

j= i 


«+j 

M 


\-x) 


§ Wr(i + ^)-r(i + ^)’ 


(3.15) 


where the second equality follows from the definition of Pochhammer symbol (13.101) . In view of 
the Gauss’s multiplication formula Bn formula 5.5.6] 


T(nz) = (2vr)( 1 - n )/ 2 n nz - (1/2) ]] V (z + 

k =o ' n ' 

with z = i + and n = M, we could further simplify (13.151) to get 

pi> _ ^ \^ A\ r(« + 1 + kM)(—xM M y 

k ~ ( ^\i) T(a + l + iM)M kM 


(3.16) 


(3.17) 


Combining (13.171) with (12.71) . it is readily seen that 




(3.18) 


which is the first identity in (|3.12l) . Note that both q^’ M ^ and Pf are monic polynomials of 
degree k. 

To show the second identity in (13.121) . we obtain from (13.71) and (13.131) that 


2nik\Uf =i r(k + ^ 


pc- 1- 

Jc-i 

M 1 


c+ioo 


M 


r(a) 

F(*-fc)tV 




M \M 


-Ms 


ds 


J/[{k+l)M 


27rifc!r(a + 1 + kM) J c _ io0 T(s — fc) 


l 


c+ioo 


P(s) 


-r(Afs + 1 — M + a)x 


—Ms 


ds, 


where we have made use of (13.161) again in the second step. This, together with (12 .14 j) . implies 
that 


n v ( X M 1 

l M M J M{k+i) M ~ 


= x e 


a— x(.ct,M) 




as desired. 

Finally, the relation (13.14[) follows immediately from a combination of (11.41) . (13.51) and (|3.12l) . 
Alternatively, this relation can also be checked directly from the double contour integral repre¬ 
sentations (11.111) and (13.111) with the help of multiplication formula (|3.16l) . 

This completes the proof of Theorem 13.11 □ 


Remark 3.1. By setting x = y in (|3.14|) . we simply have that K „is related to Kf via an 
M-th root transformation. Let n —>• oo, this in turn provides the other perspective to explain 
the appearance of Fuss-Catalan distribution in biorthogonal Laguerre ensembles, since it is 
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well-known that the Fuss-Catalan distribution characterizes the limiting mean distribution for 
squared singular values of products of random matrices [3 HUD]- As a concrete example, we may 
focus on the case 9 = M = 2. According to [551 .S3], the empirical measure for scaled squared 
singular values for the products of two Ginibre matrices converges weakly and in moments to a 
probability measure over the real axis with density given by 


V3 


4 2 

2 3 7TX 3 



x e 


27 

°'T 


(3.19) 


On the other hand, by [36] (see also USD, the limiting mean distribution for scaled particles 
from biorthogonal Laguerre ensembles (11.11) with 9 = 2 takes the density function given by 


\/3 

i 

2ttx 3 



x € ^0,3^ . 


(3.20) 


Clearly, the density (13.201) can be reduced to (13.191) via a change of variable x —>• 2y/x, as 
expected. 

Remark 3.2. From [331, Theorem 5.3], it follows that, with K" defined in (|3.5I) and v\,..., um 
being fixed, 

lim (3.21) 

ti^-og n \n nJ 

uniformly for x, y in compact subsets of the positive real axis, where 


K u (x,y) 

wL 


(2m) 

r*i 


— 1 / 2+200 r 

ds / 

— 1/2—zoo */S 


M 


chf] 


,t ni —s —1 


T{s + Vj + 1) SU17TS x t y 
T(t + Vj + 1) sin7rt s — t 

3=0 j 


- f r 1,0 

Jo 


-VQ, —Vi, . . . , -I'M 


ux G. 


M, 0 
0 ,M +1 


Vi, ■ ■ ■ , VM, VQ 


uy d u, (3.22) 


and where X is a contour starting from +oo in the upper half plane and returning to +oo in 
the lower half plane which encircles the positive real axis and Ret > —1/2 for t € X. This fact, 
together with our relation (13.141) and the hard edge scaling limits of Borodin (11.51) , implies that 


Mx a I Ja +1 i (ux) Jg+i m(( uv) AI )u a du = 


x 


M -1 


M M ~ 1 


K l 


f y M x M \ 


(3.23) 


The identity (13.231) was first proved in [323 , where the authors gave a direct proof by noting that 
Wright generalized Bessel functions J a ^ defined in ([1.61) can be expressed in Meijer G-functions 
if b is a rational number. Since it is easily seen from (13.221) and the multiplication formula (13.161) 


that 


x M 1 0 ( y M x M 

M M ~ l [W 1, At™ 


K^ M \x,y), 


(3.24) 


where K^ a,M \x,y) is given in (II. 14|) . the proof presented in [32] also gives a direct proof of 
identity (11.181) if 6 = M € N. To show (11.181) for general 9 > 1, we first observe from (11.161) . 
the residue theorem and m that 


= X 

k =0 


(_l)fc Qx ke 
kl T(a + 1 + k6) 


9J a+ i,e{x 9 )- 


(3.25) 
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Similarly, by deforming the vertical line in (11.151) to be a loop starting from —oo in the lower 
half plane and returning to —oo in the upper half plane which encircles the negative real axis, 
we again obtain from the residue theorem that 



(3.26) 


A combination of the above two formulas, (11.141) and (11.171) gives us (11.181) . 


A The Meijer G-function 

For convenience of the readers, we give a brief introduction to the Meijer G-function in this 
appendix, which includes its definition and some properties used in this paper. 

By definition, the Meijer G-function is given by the following contour integral in the complex 
plane: 



where T denotes the usual gamma function and the branch cut of z “is taken along the negative 
real axis. It is also assumed that 


• 0 < m < q and 0 < n < p, where m, n,p and q are integer numbers; 

• The real or complex parameters a\,... ,a p and b\,... ,b q satisfy the conditions 


&k ~ bj 1,2, 3,... , for k = 1, 2,..., n and j = 1,2,..., m, 


i.e., none of the poles of T(bj +u), j = 1,2,... ,m coincides with any poles of T(l — a*,— u), 
k = 1, 2,..., n. 

The contour 7 is chosen in such a way that all the poles of T{bj + u), j = 1 ,... ,m are on the 
left of the path, while all the poles of T(1 — — u), k = 1 ,... ,n are on the right, which is 

usually taken to go from — ioo to ioo. For more details, we refer to the references 133 E]. 

Most of the known special functions can be viewed as special cases of the Meijer G-functions. 
For instance, with the generalized hypergeometric function p F q given in (13.91) . one has m 
formula 16.18.1] 


<2 



(A.2) 


k =1 


This, together with the fact that 



(A.3) 


gives us 



7 


(A.4) 
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